mdesc <- function(x, mean = FALSE, digits = 2L) {
  if (isFALSE(mean)) {
    sum(is.na(x))
  } else {
    round(mean(is.na(x)), digits = digits)
  }
}